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Direct photon production in Au+Au collisions at y / ijvjv=200 GeV has been measured in the 
Pt range of l<pT<5GeV/c. The yield as a function of centrality is in agreement with 
published data of the RHIC 2002 run, and with results from a new method that explores 
very low mass dileptons. The result is compared to several theoretical calculations, and it 
is found that the measurement is not inconsistent with calculations including thermal photon 
contributions. 



PACS: 24.85.+p, 25.75.-q, 25.75.Nq, 25.75.Dw 

1 Introduction 

Several discoveries from the first five years at the Relativistic Heavy Ion Collider (RHIC) indicate 
that a new hot and dense partonic matter have most likely been produced in central Au+Au 
collisions at y / s]v~/v=200 GeV [1^1]. The drastically increased hard scattering cross-section at 
the RHIC energy makes it possible to use hard probes (jets, heavy quarks, etc.) as "measures" of 
the colliding system. High px hadrons, i.e. fragments of jets, have been found to be suppressed 
as a consequence of a final state interaction (energy loss) of hard scattered partons with the 
partonic medium. This observation is supported by the high pt direct photon measurement by 
the PHENIX experiment [5]; the high p? photons originating from an initial hard scattering 
were not suppressed. The energy loss phenomena has also been observed with electrons from 
semi-leptonic decay of charm and bottom quarks [6]. In addition to new probes, conventional 
observables also provided important information. For instance, a large collective flow of bulk 
particles gives a hint of a rapid thermalization of the partonic matter, because a large asymmetry 
of pressure gradient that results in a large flow can only be realized in the early time. In summary, 
the results suggest that the matter created in heavy ion collisions at RHIC is strongly interacting 
and quite opaque like a perfect liquid. 

Provided that the hot and dense partonic matter is created, its thermodynamical nature such 
as temperature, phase transition order parameters and the degree of freedom in the medium are 
of great interest. The temperature obtained from a statistical model fit to particle yield ratios 
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is the chemical freeze-out temperature, and no longer reflecting the one in the partonic matter. 
Therefore, more direct probes are deserved. 

Photons are an excellent probe for extracting such the thermodynamic al information, because 
once produced, they do not interact strongly with a medium. They are produced through Comp- 
ton scattering of quarks and gluons and annihilation of quarks and anti-quarks (leading order), 
and bremsstrahlung or fragment (next leading order). There is also a prediction of a jet-photon 
conversion process, which occurs if the hot and dense matter is formed, by a secondary interac- 
tion of a hard scattered parton with thermal partons in the medium [7]. The contributions from 
various processes are mixed up in any px range, but each can dominantly be seen at certain £>t's. 
The calculation shown in Fig^predicts that photons with pt<1 GeV/c are mostly contributed 
by hadron gas interaction via processes of: irir(p) — > jp(w), irK* — > Kj and etc.. [7] The 
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Fig. 1. Theoretical expectation on hadron gas in- 
teraction, thermal radiation, and pQCD photon for 
central Au+Au collisions at y / sjvjv=200GeV. 



Fig. 2. Theoretical expectation including jet- 
photon conversion compared with published data 
for central Au+Au collisions at ^Js jv jv =200Ge V. 



thermal radiation from QGP state is predominant in l<pT<2.5GeV/c. For pt>2.5 GeV/c, a 
contribution from jet-photon conversion process will be seen on top of the one from the initial 
hard scattering as shown in Fig. 13 [8]. 

The PHENIX experiment at RHIC has succeeded to scope 1<pt<5 GeV/c in Au+Au colli- 
sions at v / JJv r /v=200 GeV where the thermal radiation and the jet-photon conversion process are 
likely to manifest. This paper will present a detailed description of the analysis and the results, 
together with theoretical models. 



2 Detector and Analysis 

Fig-E]shows the PHENIX detector at RHIC. Both central arms (East and West) of the detector 
include Drift Chamber (DC), Pad Chamber (PC), Ring Imaging Cerenkov Counter (RICH), lead- 
scintillator sandwich type (PbSc) and lead-glass type (PbGl) calorimeter. The East arm has a 
Time Of Flight (TOF) detector and a Time Expansion Chamber (TEC) in addition. Each covers a 
pseudo-rapidity range of \r)\ <0.35 and a quarter azimuth. A detailed description of the detector 
is given in the literature [9]. 
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West East 

Fig. 3. The PHENIX detector. 

In this analysis, PbSc was used for measuring energies of photons and 7r°'s. The analysis 
of direct photons require a precise determination of background photons decaying from known 
hadronic sources such as 7r° and 77. PHENIX has measured the transverse momentum spectra of 
7T° up to 20 GeV/c in Au+Au collisions at y / ijv r /v=200 Ge V with ~900 M minimum bias events 
from RHIC 2004 run [10]. This work used ^93 M events out of the above in order to control 
run-by-run systematic fluctuation. The momentum spectra of 77 and other hadronic sources were 
estimated by replacing px by (pt 2 — M T 2 + M?, 2 ) 1 / 2 in fit functions to tt° spectra. Normal- 
ization factors for the estimated function were determined so that 7y/7r pT=oo =0.45±0.05 [5], 
7 // 7r °p r =oo = l-0 an d w '/ 7r °j9 T =oo = l-0- r '"h e ra ti° of background photons are shown in Fig._4ja). 
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Fig. 4. (a) Ratios of 7's from hadrons to 7's from 7r°, and (b) Quantities used for photon identification. 

Photon candidates were selected by applying a threshold on a PID likelihood variable com- 
puted from several quantities such as cluster shower shape or ratio of energies among towers 
(Fig-_4jb)). it resulted in a significant reduction of hadronic clusters in the samples. The remain- 
ing hadron contamination and PID efficiency were corrected for based on a simulation study. 
Since the ratio of 7 to tt° cancels their common systematic errors, the excess of the measured pho- 
ton over the estimated background photon is evaluated in terms of (7/ ir° ) m easured / (7/ ^°)background 
(double ratio). The systematic errors are summarized in Tab.[l] The total p^-correlated error is 
7.5 %, and the total point-by-point systematic error is 7.0 %, respectively. 

In addition to the above conventional approach, a new method that utilizes very low mass 
dileptons has been applied [11]. It is well known that the e + e~ invariant mass distributions from 
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Tab. 1. Systematic Errors on direct photon and very low mass dilepton measurement 
Direct photon measurement 



7r u peak extraction 


6% 


Energy scale in terms of -f/ir° 


3% 


j/tt° background 


3% 


7 unfolding 


3% 


Acceptance in terms of j/tt 


2% 


PID error in terms of 7/77° 


1.5% 


Off- vertex contribution to 7/77° 


2.1% 


Hadron contamination estimate 


5% 


7/7T conversion correction 


1.5% 


i]/ir° estimate error 


2% 



Total pt -correlated error 


7.5% 


Total point-by-point error 


7.0% 


Very low mass dilepton measurement 


Numbers are relative error on J direct/ "fail 


h/it {) estimate 


20% 


Acceptance 


5% 


Total error 


21% 


(point-by-point only) 





7T° or r\ Dalitz decay follow the Kroll-Wada formula as shown below [12]. 

1 dN ee _2a f i^f n 2ml 1 \ F(rr) 2 ^2 n rn 2 e 3 

where m e is the mass of electrons, m ee is the invariant mass of a produced dilepton, M and 
\F( m ee)\ aie tne mass an d the form factor of a parent particle, respectively. The formula is 
also valid for internal conversion of direct photons by substituting \F(m1 e )\ with 1. The e + e~ 
invariant mass distribution of 7r° and rj Dalitz decays as well as direct photon internal conversion 
are demonstrated in Fig.|5] Using the formula, the amount of dileptons can be estimated at a given 
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Fig. 5. Invariant mass spectra by the Dalitz decay of n , rj and internal conversion of direct photons. 

invariant mass region where a ratio of the measured to the background is taken as a function 
of pr (in this analysis, 90<Af ee <300MeV/c 2 ). Since the yield at the lowest mass region is 
approximately proportional to the total number of produced e + e~, the ratio at the measured 
mass region is then converted into the one at the very low mass region (0<M ee <30MeV/c 2 ) 
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using the Kroll-Wada formula. This work used ^900 M minimum bias events. The systematic 
errors for this method are also summarized in Tab.Q 



3 Results and discussion 



Fig.|6]shows the -f/ir° double ratio for 0-10 % central collisions as a function of px- The results 
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Fig. 6. (j/ir^easured)/ (l / ^background) tot real photon measurement for 0-10 % (left), 30-60 % (middle) 
and 70-80% (right) centrality (black solid points), together with J direct flail + 1 from very low-mass 
dilepton analysis (blue open points), and NLO pQCD expectation scaled by number of binary nucleon- 
nucleon collisions (lines). 



from very low mass dilepton analysis are also plotted on the figures in blue, in the form of 
IdirectHaii + 1' which is not equivalent to 7/71"° double ratio. Both direct photons and very low 
mass dileptons are consistent given the large error of the direct photon measurement. They are 
also consistent with the published results from RHIC 2002 run [5]. In l<pT<2.5GeV/c, there 
is no significant excess in the direct photon measurement in all the centralities, while the very 
low mass dilepton has a significant excess except for peripheral events. Overlaid on the plots 
are expectations from a NLO pQCD calculation scaled by the number of binary nucleon-nucleon 
collisions [13]. The direct photon results can be well described by the expectations except the 
0-10 % centrality data, where the onset of an additional contribution is seen in A<px<5 GeV/c. 

The result for most central events is also compared with a model including the thermal radia- 
tion, and shown in Fig.0 The model assumes an initial temperature of T§ ve = 360 MeV (T™" 1 
= 570 MeV) and a formation time of tq= 0.15 fm/c [14]. The result is not inconsistent with the 
model, suggesting that there is a possibility of the existence of a thermal source. However, fur- 
ther works on reducing systematic errors on direct photon measurement, and confirmation of no 
excess in very low mass dilepton in p+p measurement are desired to make a concrete statement. 
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Fig. 7. 0-10 % central collisions compared with a model including thermal radiation from QGP. 



4 Conclusion 

Direct photon production in Au+Au collisions at v / s]vjv=200 GeV has been measured in the pt 
range of 1<pt<5 GeV/c. The yield as a function of centrality is in agreement with published 
data of the RHIC 2002 run, and with results from a new method that explores very low mass 
dileptons. The result is compared to several theoretical calculations, and it is found that the 
measurement is not inconsistent with calculations including thermal photon contributions. Fur- 
ther works on reducing systematic errors on direct photon measurement, and confirmation of no 
excess in very low mass dilepton in p+p measurement are desired to make a concrete statement. 
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